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Abstract 

We calculate the supersymmetric 0(a s ) QCD corrections to the decays % — * q 1 xf 
= 1,2) and q~i — ► g^fc = 1, • • • 4) within the Minimal Supersymmetric Standard 
Model. In particular we consider the decays of squarks of the third generation, tj and 6, 
(i = 1,2), where the left-right mixing must be taken into account. The corrections turn 
out to be of about 10%, except for higgsino-like charginos or neutralinos, where they can 
go up to 40%. 



1 Introduction 



In supersymmetry (SUSY) [p]] every quark has two scalar partners, the squarks q~i and q~R. Quite 
generally, q~i and q~R mix giving the mass eigenstates q~i and q\ (with < m^ 2 ). Whereas the 
mixing can be neglected for the partners of the light quarks, it can be important for the third 
generation due to the Yukawa coupling which is proportional to m t or m& 0. 

In the case that the gluino is heavier than the squarks, the squarks have the decay modes 
<?i,2 - > q' xf (j = 1?2), and q~i :2 — > qx® (k = 1, . . . 4), where xf and Xk denote the charginos 
and neutralinos, respectively. These decays were discussed so far only on the basis of tree- 
level calculations, for the partners of the light quarks in ref. and with inclusion of the third 
generation (tj, bi, i = 1, 2) in ref. [|J. A detailed study of the U and bi decays can be found in 

Here we calculate the 0(a s ) QCD corrections (including the exchange of SUSY particles) to 
these decays within the Minimal Supersymmetric Standard Model (MSSM). In particular, we 
do not neglect the masses of quarks and take into account q~L-qR mixing, so that the formulae 
are also applicable to the decays tj — » bx~j, U —* an d b\ — > txj, b% — > bxl- We work in 
the on-shell renormalization scheme. For the renormalization of the squark mixing angle we 
use the scheme introduced in 0, where we applied it to the case of e + e~ — ► c^g*. We also give 
a numerical analysis of the QCD corrections to the decays U — > bxj and ti — > t x^i- 

For the decays of a squark into a light quark (m q = 0) and a photino the SUSY-QCD 
corrections were already calculated in ||. The QCD corrections to the decay t — > U xl have 
been computed very recently within the MSSM in ref. [01. The SUSY-QCD corrections to the 



strong decays q — ► qg were calculated in ref. |T0 . 

2 Tree level formulae 

The squarks q~L and (Jr are related to the mass eigenstates q\ and q 2 by: 
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Their interaction with charginos xf (j = 1> 2) and neutralinos xt (k = 1, . . A) is given by 0: 

C = gq (4 P fi + ftgfc Pl) ft + <7 * (4 P « + 4 p i) *i 6i + 9 b (4 P « + 4 ^ + h - c - ( 2 ) 

q, t, b, Xj i an d Xk denote the four-component spinors of the corresponding particles. The 
respective decay widths at tree level are 

g 2 \^{m 2 q ,m 2 l , m 2 ±) 

r°(ft - q'xf) = w , Xj (t(4) 2 + (4) 2 ] X -iilklm^f) (3) 

and | 

q 2 \z (m~ , m 2 m1 ) . 

T°(ft - <ZX° fe ) = 16 ; m s Xfc ([(<4) 2 + (bl) 2 } X' - 4a»wJ (4) 

where X(x, y, z) = x 2 + y 2 + z 2 — 2xy — — 2yz and 

X = m- qi - m q , - m-±, X = m q . - m q - m-o. (5) 



The qi-q'-xf couplings and kfj are 



/v in jni "'ij *^j'2 

with 



9 - nl o%, k% = nl o] 2 (6) 



The qi-q-xl couplings a\ k and b\ k are given by 



a ik — T^in -Aim b\ k — 1Z q in B\ n (8) 

with 

A=[(?\ Bl=( h ^\ (9) 



h q r k If 

\ n RkJ 



and 



Nk2 



h{ k = Y t (N k3 sin/3 - N M cos (3) , f{ k = -^f sin 9 w N kl - (± - § sin 2 „- 
= Y t ( N k3 sin /? - N kA cos /?) , f Rk = - sin w (tan %iV fe2 - iV*i), 



Vfc2 



/i fc Lfc = -y 6 (iV fc3 cos p + iV fc4 sin /?) , f b Lk = ^§ sin %jV fcl + V2 (± - § sin 2 9 W )^^. 
h Rk = ~ Y b (N k 3 cos /? + iY fe4 sin /?) , f b Rh = ^§ sin W (tan 9 w N k2 - N kl ) . 



(10) 
(11) 
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Nij is the 4x4 unitary matrix diagonalizing the neutral gaugino-higgsino mass matrix in 
the basis 7, Z°, - H$fwP, + H$cosf3 |§ |TT|. J7 - and ^ are the 2 x 2 

unitary matrices diagonalizing the charged gaugino-higgsino mass matrix [l], [l!|. Assuming 



CP conservation, we choose a phase convention in which N^, Uij, and are real. Yf denotes 
the Yukawa coupling, 



Y t — mt/(v 2 mw sin /?), Vj, = mfe/(v2myi/ cos/3). (12) 

3 SUSY-QCD corrections 

The 0(a s ) SUSY-QCD corrected decay width can be decomposed in the following way: 

r = r° + sr® + srw + sr® + 6r g>real . (13) 

The superscript v denotes the vertex correction (Figs, lb, c), w the wave function correction 
(Figs, ld-h), and c the shift from the bare to the on-shell couplings. 5T g ^ ea i is the correction 
due to real gluon emission (Figs, li, lj) which has to be included in order to achieve infrared 
fmiteness. According to eq. (f|) ST^ (a = v, w, c) can be written as 

</ 2 A5(m|,mJ,m|±) 



Sr^q^q'xf) = * 3 ^ (14) 



[2 4 5^ + 2 k% 5k^) X - 4 m q ,m 2f (4 5k^ + k% bi% a) ) 

with X as defined in eq. @. An analogous expression holds for 8T^ a '{qi — > qxt) by replacing 
Xf - xl 4 - q,X - X', 4 - 4, k% -> 6* , «4< a > - ^ and 54 (a) - c^ a) . 
5£fj a \ Skfj a \ etc. get contributions from gluon exchange, gluino exchange, and the four-squark 
interaction. As we will see, in the renormalization scheme used the four-squark interaction 
contribution turns out to be zero. 

In what follows we will take the squark decays into charginos as example. For the decays 
Qi <? Xk one h &s to make the replacements just mentioned before. 
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3.1 Vertex corrections 



The gluonic vertex correction (Fig. lb) yields 

Si t(v, 9 ) = a. | [(4m 2 + 2X) {Cq + Ci + c 2) + (2m 2 ± +x)Cl + Bq] £ | 



37r L-- 9- ' ---- ' - ' - ' v Xf ' ' ' ' - " 

+[2m,/m s ±C 2 ]4}, (15) 



8k 



3tt 



{[(4mJ + 2X)(C + Ci + C 2 ) + (2m|± + X)d + S ] 4 

+[2m q ,m ± ±C 2 ]l%}. (16) 



So, Co, Cij an d C2 are the standard two- and three-point functions [[L3|]. In this case, B = 
Bo(m'i±,ml.,mg,) and C m = C m (m] j .,m 2 i±,mlr, \ 2 ,mg.,mi/), where we follow the conventions 

Xj " z " " z Xj " H% H 

of fiH}] . As usually, we introduce a gluon mass A for the regularization of infrared divergencies. 



The contribution to the vertex correction due to the graph shown in Fig. lc with a gluino 
and a squark q' n {n = 1, 2) in the loop is: 

(«,») = 2^ j m ^ ± [( mg , aLi? + m^cu^ - m~ g a LL ) £ q ^ + m^±a RL kij]d 

+m q >[(m q >a RL - m q a LR + m g a RR )k q ^ - m^±a LR i q ^(C x + C 2 ) (17) 
+m- g [(m q/ a RR - m g a LL + m- g a RL )k q l j - m~±a LL £ q nj ]C + (Xd + B ) a RL k q nj }, 

5k q ^ v,9) = ^^{m^±[(m q >a RL + m q a LR - m g a RR )k q nj + m^±a L R£ 9 nj ]d 

+m q/ [{m q ,a LR - m q a RL + m- g a LL ) £^ - m-^ajnk^} (d + C 2 ) (18) 
+m g [(m q >a LL - m q a RR + m g a LR ) £ q nj - m^±a RR k q nj ]C + (Xd + B ) a LR £ q nj } 



with 



CULL - [otLLjin ~ ^il ^nl> ~ ( a LR.)i n ~ T^il ^n2' ^g) 



Here, B = B (^, m^n® , m d C ro = C„(m|, mL, mj; mj, roj, ro|. 
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3.2 Wave— function correction 



The wave-function correction is given by 



5k 



~[5Zp+SZ H } & + 5Z, 



y 

q(w) 



-[SZf + SZul k% + 5Z lV k\ tj . 



(20) 
(21) 



Z , ,R are the quark wave-function renormalization constants due to gluon exchange (Fig. Id), 



bzp 9) = sz« i{9) 



2a, 

3 7T 



B + B 1 -2m 2 JB -B 1 )-r/2 



(22) 



with B m = B m (m 2 ,, A 2 ,m 2 ,), B m = B m (m 2 ,, A 2 ,m 2 ,), and due to gluino exchange (Fig. le), 
6z l1® = 2a. | cos 2^ s i + sin 2^ /jB 2 +m *[Bl + + ^isin2^(B 1 - B 2 )}}, (23) 



5Z?® - 2 -* 



{ sin 2 9q'B\ + cos 2 ? fl? + m 2 , [£} + B\ + ^ sin 2^ (£* - £ 2 )] }, (24) 

where = B m (m 2 ,, m~, m 2 -,) and = B m (m 2 ,, m~, m 2 -,). The parameter r in eq. (p2^) and 
eq. (p2|) exhibits the dependence on the regularization. As r does not cancel in the final result 
we have to use the dimensional reduction scheme [|T5|] (r = 0) which preserves supersymmetry 
at least at one-loop order. 

The squark wave-function renormalization constants Z in stem from gluon, gluino, and squark 
loops according to Figs. If - lh . They are given by: 



^(5,5) -o (^(9,9), 2 \\ rv(M) Re { S S'^( m |)} . , ., 



1 1 

m~ — m~ 

1i lil 



(25) 



with Ejj(m 2 ) = dHa(p 2 ) / dp 2 \ p 2 =m 2 . The squark self-energy contribution due to gluon exchange 
(Fig. If) is 



'3tt 



;{3B (m|, 0, m|) + 2^^?, 0, m|) 
+4m? B (m| lS A 2 , m|) + 2m| J B 1 (m|, 0, m? )}, 

and that due to gluino exchange (Fig. lg) is 



2«, 

3 7T 



D/2 2 2\ , / 2 2 2\n / 2 2 2\ 

B (m^.,m 5 , m q ) + (m & -m ? - m- g )B {m^m g , m q ) 
-2m q m g (-iy sin 26qB (m 2 q - mi m 2 ) 



(26) 
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- — m~ g m q cos29aBo[m^_, m^, m„ 



3 7T 



(27) 



The four-squark interaction (Fig. lh) gives 



«6 

67T 



sin 46*5 



A) (m|,) -AoK-J 



(28) 



where Ao(p 2 ) is the standard one-point function in the convention of JT^- Note that £®(p 2 ) 
is independent of p 2 and hence E|y (rn~ ) = £®(m~ ) = S®. 



3.3 Renormalization of the bare couplings 



In order to make the shift from the bare to the on-shell couplings it is necessary to renormalize 
the quark mass as well as the squark mixing angle: 

< (c) = sL ok se q + -R& so] 2 , 6k$ c) = sio% se, + ni so]' 2 , (29) 

54^ = Sl A q kn 56 q + 11% 5h q Rk , 5bf k (c) = Sf n B q kn 66 q + 74 bh\ kl (30) 



where 



*». = '*« = [-21 )"«■ 



(31) 



5(9*2 = Vj2/{V2nT<w sin/3) 5m t , 50^2 = Uj2/(V2mwcos{3) 5m&, and analogously for 57i| fc and 
according to eqs. (10) - (12). The gluon contribution to 5m q is 



2 a, 



<j m fo) = _J^ m [ Bo _ Sl _ r / 2 ] , S m = S m (m^, 0, 

O 7T 



2 n 



(32) 



and the gluino contribution is 



8mf 



'3tt 



m q {B\ + B 2 ) + m 5 sin 20^ - B 2 ) 



i rn 



B' 



B m (m 2 q ,m g ,m 2 g. 



(33) 



For the renormalization of the squark mixing angle we use the scheme introduced in 0: 

( g) at 8 sin A6 a 



/-> 2 2 

b7r m| x - m| 2 



Ao(m 2 )-i (m?) 



50 



(a) 



a. 



m~ g m q 



?><Kll L (m\ -m\) 



TD I 2 2 2\ ~ 2 2\ ~ 

B {m q2 ,m- g , m q ) v n - B (m^ , m s , m q ) v 22 



qv 9' q> 



(34) 
(35) 
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with Vih the Zq^ couplings, vn = 4(Jg L cos 2 6q — s^e q ) and V22 = 4(/q L sin 2 6 q — s^eq). Here 
Iq L is the third component of the weak isospin and e q the charge of the quark q. 

With this choice of 56 q the squark contribution to the correction is zero: 

Moreover, the off-diagonal contribution (i 7^ h) of Fig. lg vanishes in this scheme. We checked 

analytically that the resulting SUSY-QCD corrected decay width is ultraviolet finite. 



3.4 Real gluon emission 



The 0(oi s ) contribution from real gluon emission, as shown in Figs, li and lj, gives the decay 
width of q~i — > gq' xf- 

2Z[m|/ 00 + m\,I xl + (m| + m\, - m\±)hx + J + h\ } (36) 

where Z = [(£^) 2 + (£^) 2 ]X — 4 k\^m q tm^±.. The phase space integrals J, J n , I nm , and J™ 
have (ms., m q >,m~±) as arguments and are given in |i~4| |. An analogous expression holds for the 
Qi ~ *■ gi Xk decay width. 



4 Numerical results and discussion 



We now turn to the numerical analysis of the 0(a s ) SUSY-QCD corrected decay widths. As 
examples we consider the decays t\ — ► bxf and £1 — > tXi- Masses and couplings of the 
charginos and neutralinos depend on the parameters M, M', /1, and tan /3. The squark sector is 
determined by the soft-breaking parameters Mq, Mu, and Md, the trilinear couplings A t and 
Ab, and \x and tan (5, which all enter the squark mass matrices. For the following analysis we use 
physical squark masses and mixing angles as input parameters. We take m t = 180 GeV, mj = 5 
GeV, m z = 91.187 GeV, sin 2 6 W = 0.23, a w {m z ) = 1/128.87, and a s {m z ) = 0.12. Moreover, 
we use the GUT relations M' = f M tan 2 6 W ~ 0.5M and m 5 = M ~ 3M. Furthermore, 
we use a s (Q 2 ) = 47r/(b x) [1 - 2 6 x lnx/(6ga;)] with 6 = 11 - 2/3 n/, 61 = 51 - 19/3 n/, and 
x = ln(Q 2 /A 2 ) fL6| . Here, is the number of flavours. 
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For discussing the t\ — > b xt decay mode we take m -± = 75 GeV, m f ~ 2 = 300 GeV, m~ bi = 220 
GeV, mr = 230 GeV, cos^ = —0.4, and tan/3 = 2. In order to study the dependence on the 
nature of the chargino (gaugino- or higgsino-like), we choose three sets of M and /J values: 
M < (M = 66 GeV, ll = -500 GeV) , M ~ (M = 70 GeV, ll = -61 GeV), and 
M > \ll\ (M = 300 GeV, ll = -62 GeV). In Fig. 2 we show the dependence of the QCD 
corrections on the t\ mass (in % of the tree-level decay width) in the range of = 80 GeV 
(LEP2, Tevatron) to = 220 GeV (LHC, e + e~ linear collider). Here, we take cos% = 0.72 
which corresponds to Mq ~ My. If the xt * s gaugino-like (solid line) the correction is about 
10% to 20% for a light t x (m^ < 100 GeV) and decreases with the stop mass. In the case of 
M ~ \y\ (dashed line) the correction lies between ~ +10% and —10%. Due to the large top 
Yukawa coupling the biggest effect is found for a higgsino-like xt (dash-dotted line). Here, the 
correction lies in the range of -10% to -20%. 

The cos# f ~ dependence is shown in Fig. 3 where we plot the 0(a s ) corrected decay widths (solid 
lines) together with the tree-level widths (dashed lines) as a function of cos 6^ for = 150 
GeV and the other parameters as given above. The widths show the characteristic dependence 
on the stop mixing angle: If t± ~ tn (cos % ~ 0) it strongly couples to the higgsino component 
of xt > whereas for t\ ~ (cos6^ ~ ±1) it strongly couples to the gaugino component. Again, a 
striking effect can be seen for |/i| M. However, also for M -C and M ~ |//| the correction 
can be larger then 10%, especially if the tree-level decay width is very small. 

We have also studied the dependence on the masses and the mixing angle of the exchanged 
sbottom (Fig. lc). In the cases studied these effects turn out to be small. The dependence on 
the gluino mass is more complex: On the one hand, the gluino mass enters in the propagator 
of the graphs in Figs, lc, le, and lg. On the other hand, the mass and the couplings of xt 
in the final state also depend on rrig via the relation M ~ 0.3 m^. However, if one relaxes this 
relation keeping M fixed and varying trig the correction increases with the gluino mass. This 
has also been noticed in refs. || . 

For the decay of t\ into txi we take m^o = 50 GeV, m 4 ~ 2 = 400 GeV, m~ hi = 320 GeV, 
m^, = 340 GeV, cos 6^ = —0.4, and tan/3 = 2. Again, we choose scenarios with a gaugino-like 
Xi (M = 98 GeV, ll = -500 GeV) and a higgsino-like x°i (M = 250 GeV, ll = -55 GeV), as 
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well as one where M ~ \fx\ (M = 93 GeV, /x = —90 GeV). We show in Fig. 4 the dependence 
of the SUSY-QCD correction (in % of the tree-level width) on the t\ mass in the range of 
mi = 230 to 400 GeV, for cosf% = 0.72. As in the case of the t\ — > bxf decay, the decay width 
gets the biggest correction if Xi is higgsino-like: It is between -15% and -20% and decreases 
with increasing ti mass. For M <C (M ~ |/i|) the correction lies between ~ +8% and 
~ -5% (~ -3%). 

The dependence of the tree-level (dashed lines) and the QCD corrected (solid lines) decay 
widths of ti — > t Xi on the stop mixing angle is shown in Fig. 5 for = 260 GeV and the 
other parameters as given in Fig. 4. Again there is a strong dependence on the nature of the 
lightest neutralino. The corrections are large for a higgsino-like neutralino. 

Summarizing, we have computed the 0(a s ) SUSY-QCD corrections to the decays of squarks 
into charginos and neutralinos. We have concentrated on the third generation squarks, where 
the quark masses cannot be neglected, and the left-right mixing plays an essential role. The 
corrections strongly depend on the nature of the charginos or neutralinos. They are about 10%, 
but go up to 40% for higgsino-like charginos/neutralinos. 
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Figure captions 



Fig. 1: Feynman diagrams relevant for the 0(a s ) SUSY-QCD corrections to squark decays 
into charginos and neutralinos: (a) tree level, (b) gluon vertex correction, (c) gluino 
vertex correction, (d) and (e) quark wave-function renormalization, (f ) and (g) squark 
wave-function renormalization, (h) four-squark interaction, and (i) and (j) real gluon 
emission. 



Fig. 2: SUSY-QCD corrections for the decay ti — > bxt in percent of the tree-level width as 
a function of for m~± = 75 GeV, m t ~ 2 = 300 GeV, m bi = 220 GeV, m b2 = 230 
GeV, cos#£ = 0.72, cos0 b = —0.4, tan/3 = 2. Three scenarios are studied: M = 66 
GeV, fi = -500 GeV (solid line), M = 70 GeV, fi = -61 GeV (dashed line), and 
M = 300 GeV, fi = -62 GeV (dash-dotted line). 



Fig. 3: Tree-level (dashed lines) and SUSY-QCD corrected (solid lines) widths of the de- 
cay ti — > bxt m GeV as a function of cos$i for m^± = 75 GeV, = 150 GeV, 
m- t2 = 300 GeV, m bi = 220 GeV, m b2 = 230 GeV, cos0 b = -0.4, tan/5 = 2. Three 
scenarios are studied: (a) M = 66 GeV, fi = -500 GeV, (b) M = 70 GeV, /i = -61 
GeV, and (c) M = 300 GeV, // = -62 GeV. 



Fig. 4: SUSY-QCD corrections for the decay t\ — > t Xi m percent of the tree-level width as a 
function of for m^o = 50 GeV, m^ 2 = 500 GeV, m~ bi = 320 GeV, m^ 2 = 340 GeV, 
cos9f = 0.72, cos^ = —0.4, tan/3 = 2. Three scenarios are studied: M = 98 GeV, 
^ = -500 GeV (solid line), M = 93 GeV, /i = -90 GeV (dashed line), and M = 250 
GeV, fi = -55 GeV (dash-dotted line). 



Fig. 5: Tree-level (dashed lines) and SUSY-QCD corrected (solid lines) widths of the decay 
ti — > t xi in GeV as a function of cos % for m^o = 51 GeV, = 260 GeV, m t ~ 2 = 500 
GeV, m~ bi = 320 GeV, = 340 GeV, cos^ = —0.4, tan/3 = 2. Three scenarios are 
studied: '(a) M = 98 GeV, /i = -500 GeV, (b) M = 93 GeV, /i = -90 GeV, and (c) 
M = 250 GeV, fi = -55 GeV. 
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